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ABSTRACT 


The  well  known  analysis  of  the  single  lap  joint  by  Goland  and  Reissner 
provided  important  contributions  to  the  literature  on  stress  analysis  of  adhe¬ 
sive  joints  by  clarifying  not  only  the  importance  of  adhesive  peel  stresses  in 
joint  failure,  but  also  the  role  of  bending  deflections  of  the  joint  in  control¬ 
ling  the  level  of  the  stresses  in  the  adhesive  layer.  Subsequent  efforts  have 
suggested  the  need  for  corrections  to  the  Goland  and  Reissner  analysis  because 
of  what  have  been  conceived  as  deficiencies  in  the  model  used  to  describe  bend¬ 
ing  deflections  of  the  central  part  of  the  joint  where  a  classical  homogeneous 
beam  model  without  shear  or  thickness  normal  deflections  were  used.  Tv-o  nro- 
qar.*-  paper  addresses  the  issue  through  the  use  of  a  more  realistic  model  in 
which  adhesive  layer  deflections  are  allowed  to  decouple  the  two  halves  of  the 
joint  in  the  overlap  region  in  the  bending  deflection  analysis,  as  well  as  in 
the  analysis  of  adhesive  layer  stresses  where  such  a  decoupling  was  allowed  by 
Goland  and  Reissner.  It  is  found  that  many  of  the  predictions  of  the  Goland 
and  Reissner  analysis  are  recovered  in  the  limit  of  large  adherend-to-adhesive 
layer  thickness  ratios,  although  substantial  differences  from  the  Goland  and 
Reissner  analysis  can  occur  for  relatively  thin  adherends. 
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NOMENCLATURE 


By,  (=E't)  -  Stretching  stiffness,  upper,lower  adherent! 

Du,  Dl  v=  E'r’/12)  -  Bending  "  . 

E  -  Adherend  Young’s  modulus 

E'  [=  E/(l-v^)]  --  "  plane  strain  Young’s  modulus 

E^  -  Bond  Young’s  modulus 

Gb  --  "  shear 

HpHj  [Eq(37]  -  \j/{E'!c2jt/2},  factors  in  coefficients  of  solution  for 

1,^2  [Eq(38)]  -  Hj  divided  by  roots  U^/8  (J=l)  and  8p^  (j=2)  of  GR  solutions 

L  (^^+S/2)  -  half  length  of  joint  (see  Fig.  2) 

L  (subscript)  --  reference  to  lower  adherend 
Mu.Ml  -  moments  in  individual  adherends 

Mg  -  moments  in  loaded  adherend  at  ends  of  overlap 

R  -  ratio,  U/p 

Rj  (J=l,2)  -  ratios  of  roots  pj  to  U/(8)'^  (j-1)  or  (8)'^  P 

Tu,Tu  "  adherend  resultants,  ta^^ 

T  (x=0,2L)  -  resultants  applied  at  ends  of  Joint 

U  [Eq(A-8)]  -  t(T/Du)*^  ^  (12£j*^ 

U  (subscript)  -  reference  to  upper  adherend 
V  —  lateral  loads  at  ends  of  joint  (see  Fig.  2) 

Vg  -  lateral  loads  at  ends  of  overlap 

b  (subscript)  -  reference  to  bond  layer 

c  -  {/2  in  GR  notation 

h  -  height  (combined  thickness)  of  overlap  region  in  thickness  direction 

k  rEq(7)]  -  ratio  of  Mg  to  ft/2  or  of  displacement  at  ends  of  overlap  to  t/2 

k„s  --  k  expression  obtained  by  Hart  Smith  [3] 

k„  [Eq(41)]  -  k-"new",i.e.  k  obtained  from  current  analysis 

l4i,k22  [Eq(35)]  --  ratios  of  displacement  solution  coefficients  to  t/2 

J  (Fig.  2)  -  axial  length  of  overlap  region 

?o  "  "  --  axial  length  of  outer  section  of  adherends 

t  -  adherend  thickness 

%  -  bond  layer  thickness 

u  "  axial  displacement 

w  -  lateral  displacement 

Wy.Wu  --  adherend  displacements 

w  [Eq(14.3)]  -  (Wu+Wl)/2 

X  -  axial  coordinate 

z  "  thickness-wise  coordinate 


Abj  [Eq(22)] 

P  [Eq(A-6.3)] 
bp  lEq(14.1)] 
Ex 
7b 
X 

K 

M,,P2  lEq(22)] 
Pe 


Po 


P. 


-  coefficients  of  solution  of  &p 
(PePt)'^ 

-  resultant  difference  T^-Tl 

-  nominal  applied  strain,  a^/E' 

-  bond  layer  shear  strain 

-  dimensionless  overlap  length  t/i 

length  of  outer  adherend  segment  tjl 
“  exponential  coefficients  of  differential  equations  for  w  and  &j. 

-  Young’s  modulus  ratio  E'/E^ 

-  modulus  ratio  E'/G^, 

-  thickness  ration  t/t^ 


INTRODUCTION 


One  of  the  most  widely  quoted  papers  in  the  literature  on  stresses  in  adhesive  joints  is  that 
of  Goland  and  Reissner[li  (subsequently  referred  to  as  "GR")  on  single  lap  joints.  This  paper 
was  particularly  significant  in  being  the  first  effort  to  identify  the  effects  of  adherend  bending 
deflections  on  the  peel  and  shear  stresses  in  the  adhesive  layer  of  a  single  lap  joint. 

The  portion  of  the  GR  analysis  relating  to  bending  deflections  treated  the  actual  joint.  Figure 
1(A),  as  a  stepped  homogeneous  beam.  Figure  1(B),  in  which  the  height  of  the  center  section 
was  assumed  to  be  twice  the  thicknesses  of  the  adherends  (thus  ignoring  the  thickness  of  the 
bond  layer).  The  forces  due  to  tensile  end-loading,  assumed  to  lie  along  the  line  a-a’  in 
Figure  1(A),  were  found  to  produce  varying  moment  about  the  neutral  axis  at  any  point  along 
the  representative  beam  which  resulted  in  the  bending  deflections  of  interest.  The  deflection 
analysis  for  this  system  was  treated  essentially  by  "beam  column"  analysis  [2],  i.e.  as  the 
analysis  of  deflections  in  a  beam  with  combined  column  (i.e.  axial)  and  lateral  loading.  (Note 
that  in  the  GR  analysis  the  column  loading  is  tensile  rather  than  compressive.  The  latter 
represents  the  case  used  in  [2]  to  introduce  the  concept  of  Euler  column  buckling.) 


(1)  GOLAND,  M.  and  REISSNER,  E.,  Stresses  in  Cemented  Joints,  J.  Applied  Mechanica  (ASME),  v.l  1,  (1944),  p.AI7-A27 

(2)  TIMOSHENKO,  S.  and  GERE,  J.,  Theory  of  Elastic  Stability,  2nd  ed.,  ChapL  1,  McGraw  Hill  (1961) 
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CQ  FORCES  CONTRIBUTING  TO  MOtCNT  IN  OVERLAP 


Figure  1  Assumed  Model 
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A  question  arises  with  the  GR  analysis  in  the  manner  in  which  the  moment  is  represented  in 
the  center  section.  The  crux  of  the  matter  is  the  comparison  of  the  classical  moment-curvature 
relation  for  a  homogeneous  beam  and  the  corresponding  expression  for  a  layered  beam  (i.e. 
Equation  (1)  vs.  (2)  ): 


MOMENT-CURVATURE  RELATIONSHIPS 


HomogeneousBeam: 

,,  Eh^  d^w 

Mn  ~  — TT - T 

(1) 

12 

Layered  Beam 

(2) 

In  the  above  equations,  Mj^  represents  the  moment  about  the  neutral  axis,  w  is  lateral 
deflection  and  h  is  the  height  of  the  homogeneous  beam  --  ie.  approximately  twice  the 
adherend  thickness,  t,  for  a  thin  bond  layer.  With  z  as  the  thickness-wise  coordinate  having 
its  origin  at  the  mid  plane  of  the  bond  layer,  the  remaining  quantities  in  (1)  and  (2)  are  given 
by 

Af »  =  -  zadz 

My  =  m -2)<5^z  ,  (z  m)<5^dz 

Equation  (1)  is  the  classical  moment-curvature  relation  for  a  homogeneous  beam  of  height 
h  which  ignores  transverse  shear  and  thickness-wise  deformations.  Equation  (2),  on  the  other 
hand,  is  a  statement  of  equilibrium  of  the  components  indicated  at  the  right  end  of  the 
segment  in  Figure  1(C)  which  make  up  M[^.  Equation  (1)  happens  to  be  equivalent  to  (2) 
under  certain  conditions,  but  these  do  not  apply  to  the  ends  of  the  overlap  where  significant 
shear  straining  of  the  bond  layer  takes  place.  Because  it  provides  for  equilibrium,  (2)  is  a 
mandatory  relationship  which  it  is  desirable  to  preserve  in  the  formulation. 

It  is  of  interest  to  compare  (1)  and  (2)  for  the  case  of  a  homogeneous  beam.  Noting  the 
standard  expression  for  axial  stress  in  the  overlap  region  under  combined  bending  and 
stretching  (i.e.  a  linear  function  of  distance  from  the  neutral  axis): 

CT  =  L  -  {wheref=T,*Ty  \  /  =  /iVl2)  W 

^  2t  I  ^  ^ 

then  making  use  of  (1)  to  represent  in  (4)  and  introducing  the  latter  into  (3)  leads  to 
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(5.1)  LTl^Lt^IMn  (5.2)  LTu^Lf-ht^ 
\2  8  2  28  2  28 

Since  t=h/2,  Eq(5)  lead  to 

Mu^ML^L{TL-Tu)^lEt^^^  =  l-Eh^^^  (6) 

^  ^  2  ^  ^  3  dx^  12  dx^ 

agreeing  with  (2)  on  inserting  (1)  into  the  latter.  This  implies  that  Eq(S)  have  to  hold  in  order 
for  (1)  to  be  applicable.  In  the  single  lap  joint,  however,  as  noted  in  Fig  1(A),  and  Tl  are 
zero  at  the  left  end  of  the  overlap,  while  My  and  Ty  are  zero  at  the  right  end.  Furthermore,  in 
the  GR  analysis,  was  represented  at  the  left  end  of  the  overlap  where  (sec  Figure  2)  \=t^ 
in  terms  of  Tt/2  through  a  dimensionless  parameter  k: 

Mu(x  =  iQ)  =  hLf  (7) 

2 

where  k  was  a  function  of  load  in  the  range  0.26<k<l.  Table  1  gives  a  comparison  of  required 
traction  conditions  at  x=Cq  with  the  values  for  these  conditions  which  are  implied  by  (5)  when 
(7)  is  taken  into  account.  It  is  apparent  that  if  k  is  not  equal  to  unity,  the  required  end  conditions 
on  Ml,  Tu  and  Tl  cannot  be  realized  if  (5)  apply.  Moreover,  the  fact  that  M^  is  negative  at  x={^ 
while  My  is  positive  implies  the  presence  of  a  curvature  discontinuity  in  the  upper  adherend  if 
(1)  applies.  Thus  Equation  (1)  cannot,  strictly  speaking,  be  used  to  represent  over  the 
length  of  the  overlap.  Hart  Smith  first  pointed  this  out  qualitatively  in  [3],  concurring  with  the 
observation  that  me  use  of  (1)  requires  the  presence  of  nonzero  tractions  over  the  left  end  of  the 
lower  adherend.  Rather  than  attempting  to  retain  (1),  it  is  desirable  to  use  a  formulation  which 
is  consistent  with  (2)  regardless  of  whether  or  not  conditions  apply  for  which  (1)  is  valid. 


['^1  IiART-SMITH,  L.  }.,  Adhesive  Bonded  Single  Lap  Joints,  NASA  Contnctor  Report  112236  (1973) 


TABLE  1.  COMPARISON  OF  CONDITIONS  FROM  EQ(5)  WITH  TRACTIONS  AT 

_ ^ _ 

REQUIRED  END 

VALUES  IMPLIED  BY 

Mn 

-(l-k)ft/2 

-(l-k)ft/2 

Mu 

kft/2 

kft/2 

Ml 

0 

kft/2 

Tu 

f 

(5+3k)f/8 

Tl 

0 

(ll-3k)f/8 

(5.3) 
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As  suggested  above,  the  lack  of  validity  of  Equation  (1)  is  associated  mainly  with  the 
presence  of  shear  strains  in  the  bond  layer.  More  specifically,  the  equivalent  of  the  following 
equation  (see  Appendix  A,  Eq(A-6)  )  which  GR  used  to  obtain  the  gradient  of  shear  strain 
in  the  bond  layer: 

^  E'thu  2 

D 

gives  a  value  of  zero  when  (5)  are  allowed  for,  so  that  a  condition  of  constant  or  zero  shear 
strain  in  the  bond  is  implied  by  (5).  Lap  joints  typically  have  large  shear  stress  and  therefore 
strain  gradients  along  their  length  due  to  adherend  deformations,  contrary  to  the  latter  result. 
It  should  be  understood  that  the  impact  of  shear  straining  on  the  performance  of  the  joint  is 
associated  with  the  fact  that  organic  bond  materials  are  relatively  soft  compared  with  typical 
adherends.  The  extreme  version  of  this  situation  can  be  visualized  in  terms  of  a  beam  cut 
along  its  neutral  axis  like  a  leaf  spring.  Since  bending  stiffness  is  a  cubic  function  of  beam 
depth,  the  bending  stiffnesses  of  each  of  the  two  halves  would  be  one  eighth  of  that  of  the 
original  beam.  A  simple  analysis  shows  that  the  combined  bending  stiffness  of  the  two  half¬ 
beams  would  be  the  sum  of  the  two  components,  amounting  to  25%  of  that  of  the  original 
beam.  Thus,  shear  straining  along  the  neutral  axis  can  be  expected  to  decrease  the  bending 
stiffness  of  the  composite  beam  making  up  the  overlap  region  of  the  joint. 

In  general.  Equation  (1)  will  overestimate  the  bending  stiffness  of  the  overlapping  part  of  the 
joint,  preventing  correct  determination  of  bending  deflections.  Since  the  thrust  of  the  GR 
analysis  was  to  demonstrate  the  influence  of  bending  deflections  on  the  stresses  developed 
in  the  bond,  an  effort  to  correct  the  situation  appears  desirable.  It  should  be  expected  that 
differences  between  predictions  of  the  original  GR  approach  and  those  of  the  present  analysis 
which  provides  for  effects  bond  shear  straining  on  decoupling  of  the  two  adherends  will  be 
greater  in  joints  containing  relatively  thin  adherends.  This  is  due  to  the  fact  that  in  the  case 
of  thicker  adherends,  the  range  of  the  decoupling  effects  is  shorter  and  is  confined  to  a  region 
close  to  the  ends  of  the  overlap.  It  can  be  argued  from  St.  Venant’s  principle  that  Eq(l) 
should  be  valid  sufficiently  far  from  the  overlap  ends,  and  the  influence  of  the  load  diffusion 
process  which  makes  it  valid  will  be  less  when  the  effects  of  interest  here  are  short  in  range 
relative  to  the  adherend  thicknesses. 

In  addition  to  the  discussion  of  Hart  Smith  which  was  previously  cited  [3],  problems  with  the 
GR  treatment  of  the  bending  analysis  have  been  noted  by  othersf4-51.  Benson[4],  as  discussed 
by  Adams[5],  described  a  modification  of  the  analysis  which  dealt  with  the  situation.  In 
addition,  Hart-Smith[3]  did  an  extensive  amount  of  work  on  the  single  lap  joint  that  provided 
an  alternate  method  of  introducing  the  appropriate  bending  deflection  corrections  which  will 
be  discussed  subsequently. 

BENSON,  N.  K.  Influence  of  Stress  Distribution  on  Strength  of  Bonded  Joints  in  Adhesion, 
Fundamentals  and  Practice  ,  Gordon  and  Breach,  NY  (1969)  p.  191-205 

[S]  ADAMS,  R.C.  and  WAKE,  W.  C.  Structural  Adhesive  Joints  in  Engineering,  Elsevier  Applied  Science  Publishers,  NY  (1984) 
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DISPLACa<ENTS 


Figure  2  Joint  Geometry 

In  the  present  case,  the  center  section  of  the  joint  will  be  treated  essentially  as  a  pair  of 
coupled  beams  throughout  the  analysis.  Once  again  it  is  emphasized  that  GR  did  use  the[4] 
coupled-beam  approach  to  deal  with  of  peel  and  shear  stresses  in  the  adhesive  layer,  but  not 
in  their  deflection  analysis. 

THE  ANALYSIS  OF  GOLAND  AND  REISSNER 


General  Remarks 

It  is  useful  here  to  repeat  certain  results  which  were  obtained  from  the  GR  analysis.  Figure 


COORD  I  NATES  D 1 5PLACEMENTS 
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Figure  4  Antisymmetry  Conditions 


2  gives  the  basic  parameters  of  interest.  Here  "U"  and  "L"  denote  the  "upper”  and  "lower" 
adherends.  The  length  of  the  adherends  beyond  the  overlap  is  while  the  overlap  has_a 
length  {  (=2c,  in  the  GR  notation).  It  is  assumed  that  tension  loads  (i.e.  resultants)  T  are 
present  at  the  outer  ends.  In  terms  of  T  it  is  convenient  to  define  nominal  axial  stress,  a^,  and 


strain,  e^: 


T 

=  -  (9.1) 

t 


(9.2) 


£'  =  £/(! -v^)  (9.3) 
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(E,  V  =  Young’s  modulus,  Poisson  ratio  for  axial  straining.)  Along  with  T,  there  are  assumed 
to  be  transverse  shear  forces  V  which  produce  moment  free  conditions  at  the  end  points, 
provided 

V  =  af  ;  a  =  t/2L  ;  L  =  l^  +  H2  =  t^  +  c  (10) 

As  in  Figure  2,  the  directions  of  the  transverse  forces  are  downward  at  the  right  end  and 
upward  at  the  left  end  for  moment  free  conditions.  Fig  3  recalls  standard  conventions  for  the 
signs  of  moments,  shear  forces  and  tensions,  as  well  as  coordinates  and  displacements.  As 
usual,  positive  moments,  shear  forces  and  tensions  correspond  to  those  shown  for  the  right 
end  of  segment  "1",  while  at  the  left  end  of  a  given  segment  the  forces  and  moments  are 
obviously  the  reactions  to  those  applying  to  the  right  end  of  the  neighboring  segment  (to  the 
left  of  the  one  under  consideration).  Thus  in  Figure  3,  the  reactions  at  the  left  end  of  Segment 
"1"  are  shown  as  the  negatives  of  the  quantities  associated  with  the  right  end  of  Segment  "2". 


Additional  relationships  associated  with  the  antisymmetry  of  the  joint  are  illustrated  in  Figure 
4,  where  attention  is  focussed  on  points  located  at  and  x^  =  2L-Xl  ,  i.e.  points  located  at 
the  same  distance  from  each  end.  Equilibrium  of  vertical  forces  requires  that  the  shear  forces 
at  these  points  be  equal,  i.e. 

V(x^)  =  V(2L-x^)  (11) 


In  addition,  since  the  generic  expression  V=-dM/dx  holds,  the  fact  that  the  moment  is  a 
function  which  is  zero  at  x=L  and  whose  derivative  is  symmetric  about  x=L  requires  that  the 
moment  itself  be  antisymmetric,  i.  e. 

M(X[)  =  -M(2L-x^)  (12) 

so  that  as  shown  in  Figure  5,  the  shear  forces  at  the  ends  of  the  overlap  are  equal  while  the 
moments  are  equal  but  opposite  in  sign.  In  the  specific  case  of  the  ends  of  the  overlap,  GR 
uses  the  notation  M^,  and  for  the  moment  and  shear  force,  as  indicated  in  Figure  5.  In 
addition,  GR  expresses  M^,  in  terms  of  the  parameter  k  introduced  in  (7)  and  in  terms  of 
an  additional  dimensionless  parameter  k': 


;  V^^2k’m 


(13) 


where  X  =  Hit 

For  future  reference,  the  end  conditions  on  the  individual  adherends  may  be  stated 
as  in  Table  2.  The  last  four  rows  of  Table  2  arc  used  in  forming  end  conditions  for 
differential  equations  given  in  Appendix  A. 


Detailed  Relationships 

The  equations  encountered  in  the  GR  analysis  are  reviewed  in  Appendix  A.  Note  that 
these  exhibit  the  influence  of  lateral  deflection  on  the  moment  distribution  due  to  the 
shift  in  the  lever  arm,  illustrated  in  Figure  6.  Here  the  horizontal  component  of  end  load. 
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TABLE  2  CONDITIONS  ON  ADHERENDS  AT  ENDS  OF  OVERLAP 

II 

X 

X=«o+« 

Vu 

-V„  s  -2k'f  A 

0 

Mu 

s  ktf /2 

0 

Tu 

f 

0 

Vl 

0 

-V„=-2k'm 

Ml 

0 

-M„sktT/2 

Tl 

0 

f 

My  +  Ml 

Mo 

-Mo 

Tl-Tu 

-T 

f 

Mu-Ml 

Mo 

Mo 

Vu-Vl 

-Vo 

Vo 

T,  acts  through  a  load  line  which  shifts  with  respect  to  the  neutral  axis  of  the  adherend 
segments,  as  illustrated  in  Figure  6(A)  for  the  outer  adherend  and  Figure  6(B)  and  (C)  for 
the  overlapping  pan  of  the  joint.For  subsequent  reference  it  is  useful  to  recall  the  most  crucial 
results  given  for  the  GR  analysis  in  Appendix  A  at  this  point .  These  include  the  expressions  for 
the  quantity  k: 

Expressions  for  k 

k  = - — -  {RPTA-l.e)  ;  k'-^.  ^  UKU  {RPTA-l.l) 

where  7^,^  =tanh((/X^)  (RPTA-l.S)  ;  7^  =  tanh(f/A/2v'^)  (RPTA-7.9) 

For  1 

k  = - ! -  ;  k’=lkXU  {RPTA-%) 

1  -V^7^  4 
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together  with  those  for  the  shear  and  peel  stress  solutions: 


Expression  for  Shear  Stress 


=  cosh[v^PU-L)/r]  ^  3  (/?/T/i-10.1) 


sinh^^pX 


AX 


Maximum  Value  of  Xf^{x=L±tJ2  ;  ^PX/2>1  ;  0.262<i<l) 
0.631P  <  =v^piili£  <  1.414P 


{RPTA-\Q.2) 


Maximum  Value  of  Peel  Stress  (x  =  L±  J/2  ;  yXll  >  1 ) 

(RW^-15.3) 

or,  since  U<yin  practical  cases: 


^b^max 


Note  that  from  the  expressions  given  in  (A-6.3)  and  (A-7.2)  with  (A- 1.5),  the  following  relation 
between  the  maximum  peel  and  shear  stress  can  be  stated: 


^6^max 

T  ^ 


=  {iT. 


1  +3A: 


N 


(14) 


If  Ej/G|,  is  set  equal  to  the  expression  for  isotropic  materials  of  2(1 +Vj,)  where  is  the  Poisson 
ratio  for  the  bond  material,  the  max  peel-to-shear  stress  ratio  runs  from  about  2  for  the  upper 
limit  of  k  of  unity,  to  a  little  less  than  1  for  the  lower  limit  (k=0.26),  for  an  assumed  Poisson 
ratio  of  around  0.5.  The  specific  value  of  the  assumed  Poisson  ratio  is  not  important,  but  this 
estimate  emphasizes  that  two  maximum  stresses  take  on  values  which  are  proportional  to  each 
other  for  various  joint  geometries  and  for  which  the  proportionality  constant  is  in  the  range  of 
1  to  2. 
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My  =  V  X  +  Tw 


CB[)  MOMEKT  IN  OVERLAP  —  ND  DEFLECTION 

j^-^. - \i 


CC)  MOMENT  IN  OVERLAP  —  DEFLECTION  PRESENT 

4-w 


I - 1 

I  =VX  -  \f  I 


=  VX  - 


+  wT 


Figure  6  Effect  of  Deflection  on  Moment  Distribution 

A  certain  amount  of  controversy  has  arisen  over  the  results  given  in  Appendix  A  for  the  peel 
stresses  ( Eq(A-15) )  from  which  (A-15.3)  is  obtained.  Kuenzi  and  Stevens[6]  stated  a  version 
of  Eq(A-15)  (coefficients  of  the  peel  stress  equation)  with  the  signs  of  some  terms  reversed  from 
those  of  the  GR  version  (see  discussion  following  Eq(A-15)  )  without  commenting  on  the 
justification  for  differing  from  GR.  Moreover,  Ref[3-4,7]  claimed  to  have  verified  the  version 
given  in  [6],  while  Carpenter[8]  found  the  original  version  as  given  above  to  be  valid. 

Considerable  care  has  been  taken  here  to  identify  correct  boundary  conditions  at  the  ends  of  the 
overlap  in  order  to  be  sure  that  the  correct  expressions  for  peel  stress  coefficients  were  obtained. 
As  far  as  can  be  determined,  the  original  expressions  given  by  GR  are  the  correct  ones.  In  any 
event,  the  terms  in  question  which  correspond  to  the  term  U  in  the  second  factor  of  (A-15.3) 
summarized  above  do  not  usually  contribute  strongly  to  the  peel  stress  predictions. 

(61  KUENZI,  E.  and  STEVENS,  G.,  Determination  of  Mechanical  Properties  of  Adhesives  for  Use  in  the  Design  (^Bonded  Joints,  Foreit 
Products  Laboratory  Note  FPL-OI I  (1963) 

(7|  KUTSCIIA,  D.  and  UOFER,K.,FeasibUityofJoiningAdvancedCompositeFlightVehicle  Structures,  Air  Force  Materials  Laboratory  Report 
AFML-TR -68-391,  (1968)  p.  56,59 

(8]  CARPENTER,  W,  Goland  and  Reissner  were  Correct,  J.  Strain  Analysis,  v.  24,  no.  3  (1989) 
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IMPROVED  ANALYSIS 


Deformations  and  Shear  Stresses 

In  the  following  discussion,  it  will  be  convenient  to  represent  the  individual  adherend 
displacements  and  stress  resultants  in  terms  of  differences  and  sums,  since  the  bond  layer 
stresses  are  more  easily  represented  by  differences  in  the  adherend  displacements  and  resultants 
than  in  terms  of  their  individual  values.  Accordingly,  the  following  notation  is  introduced: 

5j<j:)=7y-r^  (15.1)  w(jc)=l((wy  +  w^)  (15.2) 


-^L 


(15.3) 


Note  that  the  individual  adherend  resultants  can  be  expressed  in  terms  of  T  and 

as  follows:  1  -  „  1- 

T^  =  ^[r^5j.(x)]  ;  r^  =  ^[r-5j-(A:)]  (16) 

In  addition,  note  that  the  quantity  which  appears  in  Eq(2)  as  well  as  (A-6.1)  for  the 

assumed  case  of  identical  adherends,  is  given  by 


My +iV/^  =  2D 


u- 


d}w 

dx^ 


(17) 


For  increased  accuracy  in  the  following,  the  quantity  t/2  appearing  in  Eq(2)  and  (A-2.2)  will 
be  replaced  by  (t+tj,)/2  to  provide  for  the  effect  of  bond  layer  thickness.  The  replacement  of 
Eq(l )  by  (2)  modified  in  this  way  can  then  be  accommodated  by  using  Eq(A- 1 .3,  .4)  to  represent 
My  and  Mj^  in  terms  of  displacements  (setting  Du=Dl)  i.e.: 

d^w,  t  +  t. 


d}w 


u 


dx^ 


dx^ 


■iTu-T^ 


or,  in  terms  of  the  notation  just  introduced, 

d^w 


dx^ 


(18) 


Then  equating  the  modified  form  of  (A-2.2)  with  (18)  (letting  the  neutral  axis  displacement  be 
represented  by  w)  leads  to 

D,;_- - l6..  =  7ajc-( _ t-w)T  (19) 

^  dx^  2  ^  2 


Note  that,  as  opposed  to  the  original  GR  analysis,  coupling  occurs  between  w  and  the  adherend 
resultants,  so  that  w  cannot  be  determined  independently  of  5p  in  the  updated  analysis.  A  second 
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equation  in  w  and  Sj-  is  obtained  by  multiplying  (A-6.1)  by  G,,  to  get  dXj/dx  on  the  left  which, 
according  to  (A-3.1),  is  equivalent  to  d^Ty/dx^,  or  equivalendy  (according  to  the  left  member 
of  (16) )  to  (d^&j'/dx^)/2.  Then  again  allowing  for  the  assumption  of  identical  adherends,  making 
use  of  (15.1)  and  (17)  leads  to 


2G.  1 

•  __i(  *  5r 

tu  Bn 


dx^ 


We  can  combine  Eq(19)  and  (20)  with  (A-5.1)  (which  there  was  no  need  to  modify)  and  obtain 
a  self-contained  system  of  differential  equations  for  w  and  5^.  After  transposing  and  re-applying 
previously  developed  notation,  we  get 

Displacement  in  Outer  Segment  of  Adherend  (jc<{^ 


dfWy  (/2 


- - - W,,  = - ccc 

dx^  <2  "  ,2 


(21.1) 


Displacement  in  Overlapping  Segment  of  Joint  +  {) 


- - — w  =  — (ax-. 

dx^  8t^  8/2 


1 


2  D 


(21.2) 


ResultantDifference  (x<l^<ll^  +  H) 


-2^Sr-2G,.' 


hdx'‘ 


(21.3) 


Eq(21.2,  .3)  can  be  solved  as  a  2x2  system.  Boundary  conditions  include  the  end  conditions  on 
5j-  in  terms  of  the  values  of  Ty  and  Ty  given  in  Table  2,  together  with  the  joining  conditions  on 
the  displacements  at  x={g.  The  homogeneous  part  of  this  system  can  be  written 


_//  c/2_  r  +  r^ 


- 

2r 


=  0 


(22.1) 


4*  -  =  0  <22-2) 

r 

where  U  and  P  which  are  key  parameters  in  the  GR  analysis  (EQ(A-8,  A-9)  are  defined  by 

U  =  t^f/Du  s/l2e,  ;  P  =  /pc/p,  ^here  pg=G^E'  ;  p,  =  tfft 

These  represent  exponential  decay  factors  in  the  functions  describing  bending  and  bond  shear 
stress  distribution  along  the  overlap  region  in  the  GR  analysis.  Solutions  to  (22)  can  be 
expressed  as  exponential  functions: 

When  substituted  into  (22),  these  lead  to  a  2x2  system  of  algebraic  equations  with  Aj^, 
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=  (23) 

Wji  and  p  as  unknowns: 

(2p2  -  C/2)£V;,  -  6(1  *p,)A^  =  0  (24.1) 

-2pV£'V^A  ■*•  (p^  -2p2) A;,  =  0  (24.2) 

The  following  notation  in  terms  of  U  and  P,  will  be  useful: 

/?=(//p  (25.1)  v.0'=1.2)=pj/p2  (25.2) 

l?i=)^Pi/£/  (25.3)  )?2  =  P2/{\/^p}  (25.4) 

where  py  are  the  roots  of  the  polynomial 

p^-[(8+A)^4^]pV-^/?^P^=0  (26) 

Pt  2 

obtained  from  setting  the  determinant  of  (24. 1 ,  .2)  to  zero. 


Note  that  we  can  express  R  given  in  (25. 1),  in  terms  of  U  and  P,  as 

R=  112. 


N 


<^x  h 


(27) 


from  which  we  observe  that  as  the  adherend  thickness,  t,  grows,  R  approaches  zero.  In  addition, 
Eq(25)  give  us,  for  Rj  and  R2 

r- 

p^=_Ly?p  (28.1)  ;  P2=v'8RoP  (28.2) 

which,  taking  note  of  (25.2),  leads  to 


«i  = 


(28.3) 


^2  = 


4.  (28.4) 


R  "  8 

In  terms  of  v  (suppressing  the  j  subscript),  Eq(26)  can  be  expressed  as 
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2 

v2-[8(l+— )+— lv+/?^  =  0  (29) 
4p,  2 


The  roots  of  (29)  are  of  the  form 


\.,\.  =  a±b  (30.1)  where  a  =  [4(ulp,)  +  (30.2)  ;  b^\la^-R^  (30.3) 

4  4 

It  is  useful  to  observe  here  that  the  standard  identity  between  the  constant  term  of  a  quadratic 
expression  and  the  product  of  its  roots,  when  applied  to  (29),  implies  that  V2  =  R^A'i-  When  this 
relationship  is  applied  to  (28.3,  .4),  the  following  relationship: 

'  (31) 

is  obtained.  For  most  practical  situations,  the  value  of  R  obtained  from  (27)  will  be  on  the  order 
of  1  or  less.  On  the  other  hand,  the  quantity  "a"  given  in  (30.2)  can  not  be  less  than  4,  so  that  in 
general  we  have  R^«a^.  Accordingly,  a  binomial  expansion  of  (30.3)  gives  the  approximation 

b^a-^  (32) 

2a 

from  which  (30)  lead  to 

(33.1)  ;  V2»2a  (33.2) 

or,  making  use  of  (28.3,  .4)  followed  by  (28.1,  .2): 

(34.1)  ;  «2“^  (34.2) 

i.e.for  large  t,for  w/iicft(p,,/?)— >0anda— >4: 

/?i»/?2=l  (34-3) 

=  ;  pj.yrp  (35) 

which  are  the  roots  of  the  differential  equations  in  the  original  GR  analysis  for  deflection  of  the 
center  part  of  the  joint  and  the  bond  shear  stresses.  In  the  updated  analysis,  the  solution  given 
in  (A-7.1)  for  the  outer  part  of  the  adherend  is  still  applicable,  while  for  the  overlapping  part  of 
the  joint,  the  complete  solution  corresponding  to  Eq(2 1 .2,  .3)  is  a  combination  of  the  exponential 
functions  appearing  in  (23)  containing  the  coefficients  ±pj  and  ±P2  (arranged  terms  of 
hyperbolic  functions  to  accommodate  the  anti-symmetry  of  the  problem)  plus  a  linear  term 
representing  the  particular  solution  of  (21.2).  In  the  following  it  is  convenient  to  express  the 
coefficients  of  terms  describing  w  as  multiples  of  t/2,  as  was  done  earlier.  Thus  the  coefficients 
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W^j  0=1.2),  where  j  corresponds  to  which  root  of  (26)  is  referenced,  can  be  conveniently 
expressed  in  terms  of  t/2  by 

(36) 

The  complete  solution  for  w  and  5^  is  then  given  by 


_  r  smh(pi(jc-L)/r] 

W  =  - 

2^^  sinh(piXy2) 


t  sinh[p2(Jc-f-)/r]  / 

+  -k^-, - -  ajr 

2  sinh(p2Ay2)  2 


(37.1) 


^  "  ^hi 


sinh[pi(j: -L)/r] 
sinh(pjX/2) 


^h2 


sinh[p2(JC-L)/r] 

sinh(p2X,/2) 


(37.2) 


It  is  convenient  to  express  the  coefficients  of  Eq(37.2)  as  well  as  (37. 1 )  in  tenns  of  k2j.  For  this 
purpose  note  that  (24.2)  together  with  (36)  imply  that  Aj,j  and  k2j  are  related  by 


A^j(J=\,2)=HjE%jl  (38.1)  where 


(38.2) 


From  Eq(9)  together  with  the  definition  of  U  following  (22),  we  can  substitute 

E't  =  12_L 
f/2 

into  (38.1)  to  eliminate  E't.  In  addition,  it  is  useful  to  express  the  H^’s  in  terms  of  the  roots 
U/(8)*^  and  (8)*^  B  of  the  GR  deflection  and  bond  sh^-ar  stress  analysis,  as  follows: 

//i  =  -7i(i^)=-7ii^  (39.1)  ;  //2=y2(8P^)  (39.2) 

o  o 

where  the  Jj’s  are  constants  to  be  determined.  Substituting  (39)  into  (38.1)  together  with  the 
expression  for  E't  just  given  then  leads  to 

(40.1)  ;  A^  =  48J2^f  (40.2) 

4 

In  addition,  substitution  of  (38.2)  into  (39),  representing  the  pj’s  in  terms  of  (28.1,  .2),  and 
allowing  for  (31),  leads  to 

1  ^2 
7i= - (41.1)  ;  J2  = - ^  (41.2) 

R?-£!  ''"a-' 

^  1  £. 
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In  the  following,  the  deflection  equation,  (A-7.1),  for  x<  will  be  expressed  in  terms  of 
revised  notation  in  which  "k„"  ( the  subscript  "n"  designating  "new")  will  be  substituted  for  the 
"k"  of  the  original  GR  analysis,  ie.  (A-7.1)  is  expressed  as 


x<t 


2  "sinhCC/X^ 


(42) 


Continuity  conditions  at  \=tg  can  now  be  established  to  obtain  equations  for  determining  Iq, 
together  with  the  k2j’s.  Note  here  that  for  practical  purposes,  the  displacements  of  the  neutral 
axis  in  the  region  of  overlap  may  be  taken  as  equivalent  to  that  of  either  adherend,  since  the 
actual  difference  is  on  the  order  of  the  adhesive  peel  strain  times  the  adhesive  layer  thickness, 
a  much  smaller  quantity  than  w^,  wl  or  w  can  be  expected  to  amount  to.  Thus  the  displacement 
of  the  outer  adherend  can  be  considered  equal  to  that  of  the  neutral  axis  at  x=lg  and  similarly 
with  the  slopes.  Continuity  requirements  then  amount  to  equating  the  value  and  derivative  of  w 
given  by  (37.1)  witii  those  of  w  given  by  (42)  for  x^C^,  [i.e.  (x-L)/t  =-X/2]  as  well  as  setting 
given  by  (37.2)  to  T  at  that  location.  After  canceling  out  common  factors  of  T  and  t/2,  three 
equations  relating  k^,  k2i  and  k22  are  obtained: 


3  J'l 
+48_ijt22  =  -1 

4  /j2  22 

(43.1) 

(43.2) 

*22  -  ‘  0 
^h22  ^hl 

(43.3) 

where  ~  tanh(pyX/2) 

(43.4) 

A  straightforward  elimination  approach  can  be  used,  here  in  which  (43. 1)  is  first  transposed  to 
get  k22  in  terms  of  k2j: 


^22  ~  ^^(^1^1  ~  ^2^ 
where 


(44.1) 


=  4&/2  (44.2) 
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which  may  be  substituted  into  (43.2)  to  eliminate  k22: 

(l+/?2Ci)yk2i-^^C2+it„  =  l 

or,  after  transposing : 


^1  “ 


1 


(45) 


Back-substituting  this  into  (44.1)  to  eliminate  k22  gives 


^2  ~  ^^(^1 


-Cf) 


(46) 


Finally,  substituting  (45)  and  (46)  into  (43.3)  provides  an  expression  for  of  the 


form 


1  +R^C2-k„ 

«,« - ^  ♦  8fij 


l+R^C^  *h22 


"  7, 

which  transposes  to  give : 


7-„,  R\C,  -Q  _ 

V6  ^  n 


'h22 


/?i(l  +/?^C2)  +  8/?2-iliil?(Ci  -C2) 


it  = 


'A22 


y?i  +  8/?2-!^^Ci  +  #"(1 

^A22 


(47) 


As  noted  previously  [  Eq(27)  J,  R  tends  toward  zero  for  large  t,  so  that  all  the  terms 
in  (47)  containing  R  as  a  factor  will  disappear,  while  R^  and  R2  approach  unity.  We 
then  have: 
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t»tb  (R,pt^O) 


1 


^*21 


'  A1 


(48) 


i.e.  k„  approaches  the  value  given  in  the  GR  analysis  for  their  "k"  given  by  (A-7.6).  (It  can 
easily  be  shown  that  of  the  present  analysis  approaches  7,^2  of  GR  analysis.)  In 
addition,  substituting  (30.2)  into  (49.2)  and  the  latter  into  (54)  gives 

t»tb  (R,pt-^0) 


A 


1 


(49) 


From  (49)  it  is  apparent  that  we  can  set 

Ji=l,  J2=1/3  (t»tb).  (50) 

and  making  use  of  the  notation  of  (44.2)  as  well  as  (45,  46)  leads  to 

t»tb  (R,pt->0) 

:  1:22  =  ^[2(i-*,)-i]  «  0  (51) 


which  is  again  consistent  with  GR.  Applying  (50,  (51))  to  (41)  then  leads  to 
t»tb  (R,pt->0) 

4*/=--2(l-U5,r  ;  A„»(|(l-iJ-l]a,t  (51) 


At  this  point  we  focus  attention  on  the  bond  layer  shear  stress  distribution.  Noting  the  relation 
between  Tjj  and  Sp  given  in  (15)  and  applying  (A-3.1)  to  (36.2)  with  expressed  in  terms 
of  &p  leads  to  the  following  for  Xj,: 

p,  cosh[ji,(x-L)/d  H2  .  coshlfijlx -£.)/<! 

*  2(  sinh(MjX/2)  2/  ^  sinh0j2^72) 

For  the  case  of  t»tjj,  substituting  (51)  into  (52)  and  allowing  for  (49)  results  in 
t»tb  (R,p ->0) 

Tj^=.rr^(  n  >11- J  coshl^P(x-L)//I  ^  3  coshlU{x-L)/^/8t]  ^ 

8  sinh(>^pX/2)  8^8  sinh(«A./2y8) 

from  which  the  maximum  value  of  shear  stress  in  the  bond  layer  is  approximately 
t»tb  (R,p,-»0) 

*  ^-2^(/(1-A,)  ’  (54) 

•h22  Sy/i  ‘h2\ 

Note  that  as  in  the  case  of  k„  given  by  (46),  Eq(54)  is  in  essential  agreement  with  Eq(A-lO) 
obtained  by  GR,  except  that  the  second  term  in  the  latter  disappears  for  sufficiently  large  \ 
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(i.e.  dimensionless  overlap  length),  whereas  here  the  second  term  persists  no  matter  how  long 
the  joint  is. 


Peel  Stresses  in  the  Modified  Analysis 

Note  that  (see  Eq(A-2.2),  Appendix  A)  the  same  expression  (  Eq(2)  )  for  the  moment 
about  the  neutral  axis,  applies  for  the  present  analysis  as  that  used  by  GR.  However,  for  the 
individual  adherends,  this  needs  to  be  examined  further  to  insure  a  correct  version  of  the 
formulation  for  the  peel  stresses.  In  particular,  the  influence  of  bond  shear  and  peel  stresses 
on  the  moments  in  the  individual  adherends  has  to  be  re-examined.  In  the  following,  it  will 
turn  out  that  the  influence  of  peel  stresses  on  the  adherend  moment  distribution  will  be 
unchanged  from  what  it  was  in  the  GR  analysis.  However,  the  contribution  of  the  shear 
stresses  on  aaherend  moments  is  influenced  by  adherend  deflections  through  the  same 
mechanism  as  that  which  applied  to  the  moments  generated  by  the  horizontal  components  of 
the  end  loads.  A  particular  point  which  should  be  understood  has  to  do  with  the  way  that 
these  horizontal  loads  generate  moments  at  various  points  in  the  adherends.  These  are  not 
large  defection  effects  corresponding  to  large  values  of  either  lateral  deflection  or  slope. 
Rather,  they  should  be  thought  of  as  long  range  effects,  since  they  correspond  to  relatively 
large  horizontal  offsets  between  the  point  where  the  load  is  generated  and  the  point  under 
consideration  where  moments  are  being  calculated.  Figure  7(A)  denotes  the  incremental  loads, 
dF^(x)  and  dFg(x)  generated  on  the  upper  adherend  by  the  shear  and  peel  stresses  at  a  generic 
axial  position  x,  where 


dF^(x)  =  ;  dFg(x)  =  Obdx 

As  stated  previously,  the  theory  developed  here  is  a  small  strain,  small  deflection  theory. 
Thus,  for  example,  rotation  will  produce  a  vertical  component  of  dF.^  given  approximately  by 
w'dF.j  at  X  which  will  produce  a  moment  contribution  at  x  given  by  (x-x)w'dF^.  The  fact  that 
w'  is  assumed  to  be  a  negligible  quantity  does  not  rule  out  the  possibility  that  (x-x)w'  may 
not  be  if  x-x  is  large  enough. 

On  the  other  hand,  a  similar  horizontal  component  of  dF^j  given  approximately  by  w'dF^j 
produces  a  moment  at  x  given  by  [w(x)-w(x)lw'dFjj;  but  since  w(x)-w(x)  is  equivalent  to 
(x-x)w')avc  where  w')ave  is  the  average  slope  in  the  interval  between  x  and  x,  the  contribution 
of  the  rotation  of  dF^  is  of  the  order  of  (x-x)w'^dFjj,  i.e.  proportional  to  the  square  of  the 
slope,  and  can  be  ignored  in  comparison  with  w'(x-x)dF^  Thus,  while  the  rotation  of  dF^,  can 
be  ignored  in  a  small-slope  analysis,  that  of  dF.^  cannot  necessarily  be  similarly  ignored. 
Accumulation  of  dF^  and  dF^  by  integration  along  the  adhesive-aaherend  interface  will  then 
prcxluce  a  moment  at  the  point  of  observation,  x,  which  is  determined  by  the  offset  of  the 
load  lines  of  incremental  loads  dF.^  and  dF^,  applied  at  x,  from  the  point  of  observation,  ie. 
by  and  shown  in  Figure  7(B).  The  adherend  moments  can  thus  be  written  as 
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Figure  7  Effects  of  Adherend  Deflections  on  Moments  in  Overlap  Region 

Mi/jc)  ^Vx*  fwfjix)  *  {xJ:)dF^(x)  -  tg{x,x)dF^(x)  ;  (x)  =  j^^^{xJ)dF^  *  J^*tg(x,i)dFg 

where  =  x  -x 

for  both  adherends,  while  expressions  for  8*;^  and  are  obtained  as  follows:  for  the  upper 
adherend,  as  illustrated  in  Figure  7(C)  and  (D),  the  load  produced  by  dF^  at  x  is  offset  from 
the  neutral  axis  at  x  by  6^  given  by 
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VpperAdherend 

=  +[wyW-H-y(i)-(j:-i)>v^i)]  (55) 

This  is  obtained  from  considering  various  contributions  to  shown  in  Figure  7.  It  is  clear, 
first  of  all,  that  a  moment  is  generated  by  dF^  from  the  difference  in  at  x  and  x,  as  shown 
in  Figure  7(C).  In  addition,  as  in  Figure  7(D)  (lower  expression  at  right  end  of  adherend 
segment),  the  tilt  of  the  tangent  vector  at  x  has  the  effect  of  counteracting  the  latter 
displacement  difference  by  an  amount  (x-x)w',  resulting  in  Eq(55)).  As  apparent  in  Figure 
7(B), however,  the  deflection  of  the  lower  adherend  tends  to  move  its  neutral  axis  toward  the 
load  line  of  dF^  rather  than  away  from  it  as  in  the  case  of  the  upper  adherend.  As  a  result, 
the  terms  in  Eq(55)  which  are  in  square  brackets  are  subtracted  for  the  lower  adherend 
rather  than  added,  resulting  in 

-  [w^(x)-w^(i)-(x-i)H'^(i)]  (56) 

insening  Eq(55)  and  (56)  into  the  integral  expressions  just  given  for  the  moments  in  the 
adherends  at  a  generic  observation  point  x,  (  x>ijj)  leads  to: 

=  V  +  fH'y(j:)  +  J^\(i)(.i(/+/^  +  [W(;(x)-w„(i)-(x-i)wy(i)]  ]<ix  *  ^  {x-i)a^{i)dx  (57.1) 

and 

-[w^(x)->v^(i)-(x-i)H'^(i))  ]dx- j^^(x-x)o^{x)dx  (57.2) 

Derivatives  of  Eq(57)  are  needed  for  the  subsequent  development.  These  are  obtained  by  use 
of  the  following  identities  which  apply  to  arbitrary  functions  F(x),  G(x)  and  H(x): 

4-["fmdx-Hix) 

dxJ*, 

from  which 

1-  rV(i)C(i)di  =  F(x)G(x)  ;  [^F(x)G{x)dx  =  G'(x)  [^F{x)dx  *G{x)F(x) 

dx-t,  dx  A,  Ji, 

leading  to 

-G{x)]F(x)dx  =  G\x)j^’‘F{x)dx  (58) 

Applying  this  to  (57)  leads  to 

=  l7+fwy(x)  +  .i(/-*-<i)Ti,(x)  -  J^*[w^x)-H'^i)lT^(i)di  -  J^*a^(i)di  (59.1) 
-H'i(i)lXj,(i)di  ♦  j^’‘a^{x)dx  (59.2) 


while  differentiating  (59)  to  get  second  derivatives  of  and  Ml  results  in 
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d^M 


dx^ 
A 


^  -  fw'^x)  *  *  w"(x)j[\(x)dx  -  a,(x)  (60.1) 


- y^L{x){^x^ix)dx  +  a^(x) 
dx^  2  •*. 


(60.2) 


Making  use  of  (16.2,  .5)  and  (30)  to  eliminate  the  integral  expressions  in  (60)  then 
leads  to 

d^M 


=  Ww"(x)  *  ‘  (/  >  ty^ix)  -  a,  (61.1) 
etc-'  2 


d^M 


^  =  T^(x)w''(x)  *  ht  *  gx'(x)  +  o,  (61.2) 
dx'^  2 


and  subtracting  (61.2)  from  (61.1)  results  in 

d\My-Mi) 


dx^ 


=  (7-yw"-rX')-2a,  (62) 


Note  the  identity 


or,  making  use  of  (16.7)  together  with  the  notation  of  (30), 


rit  rwm  //  'T'O//  O 

u^u'T^l'^l  =7^8^  + Sr'*' 


r.n 


(63) 


On  allowing  for  (14.3,  .4)  and  (17.2,  .4),  the  second  derivative  of  the  difference  between 
adherend  moments  can  be  expressed  as  : 


d\My-M0  „  A 


^D, 


dx^ 


dx* 


^Dnt 


d*e, 


U*b- 


dx‘' 


E,  dx* 


(64) 


while  from  the  same  relationships,  the  right  hand  side  of  (63)  can  be  expressed  as 


sn  .  5  -//  ^  ^b 


E,  dx^ 


+  5r’*'  (65) 


Inserting  (65)  into  (63)  and  the  result  of  that  substitution,  together  with  (64),  into  (62)  then 
results  in 


dx*  2j2  dx^ 


(66) 


The  last  equation  is  comparable  to  (20.1)  obtained  by  GR  for  the  peel  stresses,  with  the 
addition  of  the  second  order  term  on  the  left  and  the  forcing  function  on  the  right.  The 
homogeneous  solution  for  (66)  will  be  of  a  form  similar  to  (26)  for  the  GR  peel  stress 
analysis,  except  for  minor  changes  having  to  do  with  the  fact  that  the  coefficients  for  the 
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trigonometric  factors  in  (24)  will  be  slightly  different  from  those  for  the  factors  involving 
hyperbolic  functions.  In  practice,  the  changes  are  so  small  as  to  be  of  no  practical 
significance,  so  that  (26),  with  an  appropriate  modification  of  k,  can  be  used  for  the 
homogeneous  solution  for  all  intents  and  purposes. 


The  presence  of  a  non-zero  function  on  the  right  side  of  (66)  requires  the  addition  of  a 
particular  solution  to  the  homogeneous  solution.  Obviously,  the  factors  from  which  the  right 
hand  side  is  formed  are  obtained  from  the  solutions  for  w  and  given  in  (51).  The  details 
of  this  part  of  the  solution  will  be  provided  in  Appendix  B.  Again  it  will  be  found  that  in 
practical  situations,  the  differences  introduced  by  these  terms  from  those  of  the  GR  peel  stress 
expression  (except  for  the  replacement  of  k  by  k„)  are  essentially  negligible. 


COMMENTS  ON  THE  HART-SMITH  ANALYSIS 
To  date,  the  only  serious  attempt  to  correct  the  deficiencies  of  the  GR  deflection  analysis, 
other  than  the  present  one,  appears  to  have  been  that  of  Hart-Smith[3]  which  was  developed 
on  a  NASA  contract  in  the  early  70’s.  Those  results  have  been  frequently  quoted  in  reports 
and  publications  which  have  appeared  since  then,  and  a  comparison  with  the  present  results 
is  in  order.  The  most  crucial  part  of  the  Hart-Smith  analysis  is  the  expression  given  for 
adherend  bending,  which  we  may  compare  with  Eq(21.2)  of  the  present  analysis.  This  is 
repeated  here  for  comparison: 


d^w 

dx^ 


8/2 


+ 


2 


{RPT2\2) 


The  comparable  equation  in  [4]  is  equivalent  to  a  twice  differentiated  version  of 
Eq(35.2) 


d^iv 

IF 


1/2  d^w 
2t^~dF 


but  lacking  the  second  order  term  which  is  needed  to  allow  for  the  influence  of 
deflection  on  moments,  ie. 


d*w 


*/)-Lt 

Du 


The  absence  of  the  second  derivative  term  here  completely  changes  the  character  of  the 
solution.  As  a  result,  the  influence  predicted  by  this  approach  gives  much  longer  range  effects 
than  those  of  the  present  analysis.  In  effect,  this  prevents  St.  Venant’s  principle  from  being 
satisfied  in  the  sense  of  not  allowing  Eq(l)  to  be  applicable  at  any  point  in  the  system,  of 
either  GR  or  the  present  analysis  is  in  the  variation  of  k  which  is  obtained.  Table  3  compares 
the  expression  for  k  vs.  overlap  length  obtained  from  the  three  approaches. 


Considering  k  as  a  function  of  UV2  s(l2e^)*^l!/2t,  i.e.  dimensionless  overlap  length  together 
with  square  root  of  the  loading  strain,  the  GR  analysis  has  the  well  known  asymptotic  limit 
k»0.262  for  large  UX/2.  The  present  theory  will  likewise  be  found  to  produce  an  asymptotic 
limit  but  one  which  may  vary  to  some  extent  from  0.262,  especially  for  thin  adherends.  In 
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TABLE  3  COMPARISON  OF  EXPRESSIONS  FOR  k  vs.  0X72 
{  case  of  tanh(UX^)»l  } 


contrast,  the  Hart-Smith  expression  has  lower  limit  of  zero  for  increasing  UX,.  Note  that  he 
latter  result  corresponds  to  zero  shear  and  peel  stresses.  This  type  of  behavior  is  not 
consistent  with  what  is  found  in  most  problems  of  load  diffusion  which  is  what  we  are 
dealing  with  near  the  end  of  the  overlap,  and  a  lower  limit  on  the  stresses  which  persists  even 
for  indefinitely  long  overlaps  and  large  loads  seems  to  be  more  in  line  with  practical 
expectations. 

NUMERICAL  RESULTS 

The  degree  to  which  the  results  of  the  original  Goland  and  Reissner  analysis  differ  from  those 
of  the  updated  analysis  of  the  present  paper  depends  primarily  on  the  R,  the  ratio  of  U  to  p. 
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GOLAND  &  REISSNER  [11 


HART-SMITH  [3) 


Figure  10  Comparison  of  Stress  Predictions  -  Goiand  Reissner  vs.  Present  Theories 

Figure  8  gives  a  comparison  of  the  GR  and  Han-Smith  predictions  for  k  vs.  UX,  while  Rgure 
9  shows  a  similar  comparison  between  GR  and  the  results  of  the  present  analysis  for  k  vs. 
UX. 

The  effect  of  R  on  the  curve  of  k  vs.  UX/2  is  effectively  brought  out  in  these  results.  It  is 
clear  that  the  largest  difference  between  GR  and  the  present  analysis  occurs  for  small 
adherend  thicknesses. 

Figure  10  shows  the  effect  of  adherend  thickness  on  the  shear  and  peel  stress  predictions 
obtained  from  the  present  analysis  for  large  enough  UX  so  that  effects  of  joint  length  have 
settled  out,  i.e.  with  k  having  reached  its  lower  limit.  For  the  case  considered  here, 

E'  =  20x1  Og  psi,  Gb=150  ksi,  Ej,=500  ksi,  tj,  =  0.008  in,  5^  =  20  ksi. 

It  is  interesting  to  note  that  the  ratio  of  maximum  peeling  to  shear  stresses  predicted  in 
Eq(15)  for  the  GR  analysis  is  close  to  unity,  so  that  the  GR  curves  lie  nearly  on  top  of  each 
other.  Considerable  deviation  from  the  GR  curves  is  noted  in  those  of  the  present  analysis. 
In  particular,  the  shear  stresses  obtained  from  the  present  analysis  are  greater  than  those 
predicted  by  GR,  while  the  peel  stresses  are  somewhat  less  than  those  obtained  by  GR. 
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CONCLUSIONS 


In  the  context  of  the  use  of  beam  theory  for  modelling  adherend  response,  the  approach 
employed  in  this  paper  appears  to  provide  the  first  completely  consistent  model  for  the 
analysis  of  single  lap  joints.  It  is  emphasized  that  the  original  Goland-Reissner  model,  in 
analyzing  the  deflections  of  the  joint,  used  a  classical  homogeneous  beam  model  for 
calculating  deflections  which  was  not  consistent  with  the  bi-layer  model  used  by  them  to 
predict  bond  shear  and  peel  stresses.  Due  to  the  long-range  nature  of  the  effect  of  bond  shear 
strains  on  the  bending  stiffness  of  the  joint  when  the  adherends  are  thin,  the  largest 
differences  between  the  present  approach  and  the  GR  approach  occur  the  case  of  thin 
adherends.  On  the  other  hand,  the  numerical  results  which  have  been  obtained  demonstrate 
that  the  original  GR  approach  produces  numerical  results  which  do  not  substantially  differ 
from  those  of  the  present  analysis.  It  is  clear  that  for  large  adherend  thicknesses,  the  effects 
which  are  ignored  in  the  original  GR  approach  are  so  concentrated  near  the  ends  of  the 
overlap  that  their  consequences  are  not  significant.  Thus  for  thicker  adherends,  the  GR 
approach  can  be  considered  to  have  been  satisfactory  from  the  stand-point  of  sound 
engineering.  The  insight  that  Goland  and  Reissner  brought  to  bear  on  the  effects  of  adherend 
bending  deflections  together  with  the  role  of  peel  stresses  on  the  strength  of  adhesive  joints 
still  stand  as  a  major  advance  in  the  understanding  of  factors  important  to  successful  adhesive 
joint  design. 
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APPENDIX  A  EQUATIONS  OF  GOLAND  REISSNER  ANALYSIS 


BASIC  RELATIONSHIPS 
Adherend  Constitutive  Relations 

dui,  dU[  d^w.j 

dx  dx  dx^ 

(A-1.5);  (A-1.6) 

E*  -  £/(l  -v^)  ;  V  s  Poisson' sratio 

Moment  Distribution  (Figure  6) 

OuterEndoJUpperAdherend(  0<x<t^) 

Mjj^f{(xx*w)  (>4-2.1) 


M, 


d^w, 


dx^ 


(A-1.4) 


Moment  about  Neutral  Axis  in  Overlap  Region(t^<x<lg+t) 

M^  =  fax-(^-  w)f  (A  -2.2) 

Equilibrium  Relations 

Upper  Adherend  Forces  {FigA-\{A)) 

dT  dv 

(A-3.1)  i.e.  Tu^f*^\dx'  (A-3.2)  ;  (>4-3.3)  i.e.  Vy  =  (A-3.4) 

dhdjj  I  t 

(A-3.5)  ;  d^  =  ^x,-V,  (A-3.6) 

Global  Equilibrium  in  Overlap  Region  (FigA-\{B)) 

=  f  (A-3.7)  ;  (A -3.8) 

Bond  Laver  Constitutive  and  Kinematic  Relations 

Bond  Layer  Strains  (FigA-l) 

=  (A -4.1)  ;  (A -4.2) 

'b  * 

BondLayerStresses 

(A -4.3)  ;  ^b'^b^ib  (A -4.4) 

Note  that  Eq(A-2)  exhibit  the  influence  of  lateral  deflection  on  the  moment  distribution  due 
to  the  shift  the  lever  arm  through  which  the  end  load,  T,  acts  with  respect  to  the  neutral  axis 
of  the  adherend  segnments  under  consideration  (see  Figure  6(A)  for  the  outer  adherend 
and  Figure  6(B)  and  (C)  for  the  overlapping  part  of  the  joint. 


-29- 


X 


-rbCbc 


I-  q,cix 


(A)  Differential  Relationships 

'£F,^Vl*i-Vl,)-a,dx^O  ^-0 


5^  My  (about p)sM^  -My  *Vy  <fe  -j-x^ilx 


dx 


V 


u 


J2Mi^iaboutp')  =  Ml-Ml*Vj^dx-lx^dx 


dx 


(B)  Global  Relationships 


Figure  A-1  Equilibrium  Relationships 


E 


v^.v^  =  -v 
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CBJ  EFFECT  OF  ADHffiENO 
nOTATION 


Figure  A-2  Bond  Layer  Shear  Strain  Components 
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DIFFERENTIAL  EQUATIONS 

In  the  following,  the  relationships  on  which  the  diferential  equation  under  consideration  is 
based  are  given  preceeding  the  DE’s. 

Adherend  Deflections 


Outer  Segment  of  Adherend  {£^/4-1.3),(i4-2.1)) 

d^W,,  //2  f/2 

(A-5.1) 

dx^  r  r 


Overlapping  Segment  of  Joint  {£^(l),(i4-2.2)) 

(^-5.2) 

(tP’  Sp  SP  ^ 


where 


Bond  Laver  Shear  Stress  &  Strain 


U^t  _L  =/i^  (/t-5.3) 

\  Dwi 


Shear  Strain  {£^(i4-4.1),(i4-l.l)  to  {A-l  A) differentiated,  (i4- 1.5) -(A -1.6) } 

1  T’l  t  ^11 

(A -6.1) 

Shear  Stress  (£^(A-4.3),(A-1)&(A-6.1)  twice  differentiated,{A~'i.5, .6), {A-3.S) differentiated] 

_ i-SP  (A-6.2) 

dx  {2  dx 


where  P  =  Jpc/p,  I  Pg^—t  *  P/“—  ('4"6-3) 

'  p'  t 


Bond  Laver  Peel  Stress 


{Eq(A  -1 J,  .4),  (A  -3-3),  (A  -3.5,.6),  (A 


dot  -A 
__^*4Xa-  = 
dbP  P  " 


0  (A-7.1) 


£.  IM 

where  Y  =  (2 _ _)  {A-1.1) 
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SOLUTIONS 

Adherend  Deflections 


In  the  following,  the  outer  adherend  deflection  is  represented  as  a  multiple  of  kt/2, 
corresponding  to  the  notation  for  given  in  Eq(A-l).  The  parameter  k  is  determined  by 
satisfaction  of  joining  conditions  between  the  outer  adherend  and  overlapping  part  of  the  joint 
at  x=C„. 


Generic  Expressions 


x<L 


w  ■■ 


let  sinh(UxJt) 
T  sinh(t/A.J 


-cu  (/t-7.1) 


sinhl  iU/^)(x-L)/t]  _ 

sinh(£/V>/8]  2 


Continuity  Conditions  at  x=lg 

Displacement:  lcS  =  -W2*—-ai^  (A-73)  ;  Slope:  ^  ^  _  U  ^2  1  (^.74) 

2  2  ^  ^  ^  ‘  'J'h2 

Resulting  Expressions 

- Im _  (>4-7.6)  ;  l.t-W  {A-l.l) 

2  47-*/ 

where  7^,;  =tanh(t/XJ  (>4-7.8)  ;  7^  =  tanh((/X/2v^)  (>4-7.9) 

Note  that  GR  assume  the  outboard  pan  of  the  adherend  to  long  e.iuufeh  10  make  UX»1, 
leading  to  Tj,j  =1,  as  a  result  of  which  they  use  the  contracted  expressions 

it= _ ! _  ;  k'=lkUi  (>4-8) 

for  k  and  k'. 

Bond  Laver  Shear  Stress 
Generic  Solution  of  Eq(A  -6.2) 

/8  Vs"  pifiH 

Tj,=c,c  '  ♦Cjf  '  *  const.  (>4-9) 
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Bond  Laver  Shear  Stress  (con t) 

Boundary  Conditions 

(1)  Using  Eq(A  -6. 1)  multiplied  by  to  generate  dxjdx^end  conditions  determined  by  My,  My,  Ty  and  Ty, 

inserted  on  RHS  of  (A -6.1). 

(2)  Determine  const,  in  {A  -9)  by  integrating  then  inserting  into  (A  -3.2)  to  sati^  Ty=0  condition  at  x^L  +1 

Final  Form  of  Solution 

^  sinhv^px  ^ 

Maximum  Value  of  x^(,x=L±t/2  ;  y^p\/2>l  ;  0.262<JI:<1) 

0.63IP  <  a65{J!ifl}oy2piilH  <  1.414P 
a,  4 

Peel  Stresses 


(A-lO.l) 


(A -10.2) 


Functional  Notation 

Cc{x)  =  cosh(Yi_^l)  cos(7.f_il)  (A  - 1 1 . 1 )  Ss(x)  *  sinhCyiLfl)  sin(Y^l£)  (A  - 1 1 .2) 

(  t  t  t 

Scix)  =  sinh(y—)cos(y—)  (A -1 1.3)  Cs(.x)  =  cosh{y—)sin(y^)  (A -11.4) 
It  t  i 


Derivative  Relationships 

(Cc)'  =  I(5c-Cj)  (4-12.1)  (5j)'  =  I(5c+a)  (4-12.2) 

t  t 

(CcY'^-lJ^Ss  (4-12.3)  (Ss)"^2^Cc  (4-12.4) 


Generic  Solution  (ct^  .  Cf, ,  -  -constants  determined  by  boundary  conditions ) 

ai,{x)  =  aj.Cc  +  OjSs  (4-13) 
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Boundary  Conditions 


( 1 )  Eq{A  -4 .2)  combined  with  (4-4.1)  differentiated  twice  (making  use  £>/(4  - 1 . 3 , ,  4)) : 


d\  ^  E, 
dx^ 


(4-14.1) 


(2) Differentiation  o/  (4  - 1 4 . 1 ) .  making  use  o/ (4  -3. 5 ,  .6) : 


dx^  tb 


(4-14.2) 


(For  x=fg  and  {+{„,  differenced  quantities  appearing  on  the  RHS  of  Eq(A-14)  are  taken  from 
the  last  two  rows  of  Table  2.) 


ResultingExpressions 

a^(x)  =  ^III  *  ^‘P,)Cc(x)  *  (/>,^  +  |jt'/>3)5i(x)l  (4-15.1) 

where 

P\  Pi  ~  ~  ^^\x-L)lfXn  ^3  ■  ^^(x-l.)/f\a  ~  ^‘^(x-L)/fXn 


A  =  (PxP^  - P-J*i\x-L)ifxn. 


(4-15.2) 


Maximum  Value  of  Peel  Stress  (x-LtHl  ;  yX/2  >  1 ) 

2k^  kU 

Note  that  ( Eq(A  -8) )  "  =  _ 

\  2 


.■.For..L±XZ  . 


I,Tl 


1  T 


p.p.  *p.p^  p\  *pI 

AlsoJory‘Kll>\  ;  -LA _ Ll  =  _2 _ 1=1 


(4-15.3) 
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APPENDIX  B 

PEEL  STRESS  SOLUTION  IN  UPDATED  ANALYSIS 


Notation  pertinent  to  this  section  arc  given  in  Table  B.l.  Eq(66)  which  defines  the  peel  stresses 
is  repeated  here  for  reference: 


j/2  d\ 
dx^ 


+  =  Q(x) 


(5-1.1) 


where 


e(x)  =  5^(x)w"(jt)  (fl-1.2) 

It  is  noted  that  the  factors  included  in  the  function  Q(x),  on  the  right  side  of  (B-1.1),  involve 
only  the  homogeneous  parts  of  the  solutions  for  &j-  and  w(x),  so  that  (B-1.2)  can  be  written 


Q(x)  = 


sinh[^l^(x-L)/<l 

sinh(p,A/2) 


sinh[m(x-L)/r]  ^ 
sinh(}ijX/2) 


^  ,  sinh[fi,(j:-L)/d  ^  ,  sinh[pj(x-L)/d 

■7*^' 2  sinh(p,X./2  *  sinh(pjA./2) 


Defining  M  together  with  the  q^j’s  (i,j  =  1,2)  as  in  Table  B-1  and  noting  that  from  the  definition 
of  M  and  U,  we  have  that 


MU^ 


(5-3) 


then  using  appropriate  identities  for  the  hyperbolic  functions  in  (B-2)  together  with  the 
functional  notation  in  Table  B-1  leads  to 


F„(^)-F„(^) 


nf-j\  -  iT  r  ^  n  ^  -  21  12'^'  1  /n  m\ 

^  'T  *  ‘*‘^"F„(A./2)-l  ‘^”F„(V2)-1  ^■"f„(X/2)-f,j(X) 


Equation  (B-1.1)  can  then  be  expressed  as  a  sum  of  terms  for  each  of  which  the  RHS  of  (B- 
1.1)  is  an  exponential  function,  giving  a  differential  equation  of  the  type 


dx* 


dx^ 


t* 


i* 


(5-5) 


where  a  and  F  arc  constants.  Eq(B-5)  has  the  particular  solution 


a  = _ ! _ Fe"" 

'  a.* 

2?  ^ 


(5-6) 


Decomposing  the  right  side  of  (B-1.1)  into  exponential  terms  and  repeatedly  applying  (B-6) 
gives  the  particular  solution  that  is  needed  here.  The  homogeneous  solution  of  (B-1.1)  consists 
of  terms  of  the  form 


Oj  =  a*,e’’-'"cos(xx//)  +a„e’'''"sin(Y,x//) 

*o^e  ■''''cos(Y,x//)  +  <3^e  '^'"sin(  Yx/i)  {B  -7) 
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TABLE  B-1  NOTATION  FOR  PEEL  STRESS  EXPRESSIONS 


PARAMETRIC  NOTATION 
General 


T  - 


(x-L) 


Definition  of  gij’s 


.  Pi  ,  t 

2  .  A  ‘ 


M. 


<lxi' 


A  L  ‘ 


<?J2' 


A  ‘ 

8<* 


A  t 

i^. 

r  r  2 


/.e.  q,,M  -A  *  A,  k 


QijS  IN  TERMS  OF  NOTATION  OF  EQ(393,  .4),  (52.2),  (54) 


4^2%  3 


Definition  of  D,jS 

n  2f  “  2/  '*  2/  "  2/ 


FUNCTIONAL  NOTATION 

F„(^)=cosh2ji,^  ;  /='jj(^)  =cosh2pj^  ;  F,j(^)  =cosh(ji, ;  Fj,(0  =  cosh()i, +Mj)^ 
G„=sinh2>i,^  ;  G„(^)  =  sinh2m^  ;  G^^C^)  =  sinh(}i,+fi^^  ;  G„(^)  =  sinh(p, +jg^ 
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in  which  ±(Y^Yj)  [i=(-l)''^]  are  complex  roots  of  the  quartic  polynomial 

2«-t/V+4Y‘»0  (B-8) 

associated  with  (B-1.1).  The  solution  to  (B-8)  is 

z*^[±(Y,±i-7.)]*  =  ±[^±i  (B-9) 

^  \  4 

Because  of  the  second  derivative  appearing  in  (B-1.1),  it  is  found  that  y,  Yi.  contrast  with  the 
original  GR  solution.  However,  the  contribution  of  terms  involving  U  in  (B-9)  compared  with 
those  involving  y  is  miniscule,  and  can  be  ignored  with  no  risk  of  significant  inaccuracy.  Thus 
the  homogeneous  part  of  the  peel  stress  solution  is  essentially  in  agreement  with  that  of  the  GR 
solution  will  now  be 

where  the  homogeneous  term,  Cbh,  is  now  assumed  to  have  the  same  form  as  that  of  Eq(A-14) 
of  the  GR  analysis  of  Appendix  A 


but  with  Cj  and  a,  appropriately  modified: 

a,._^ - 3_f -  ;  a,._J - -  (B-U) 

in  which  Mq  and  Vg  are  as  before  given  in  (A- 13)  which  may  be  re-expressed  here  for 
convenience,  using  notation  which  has  been  introduced  in  Appendis  A  (assuming  T^,®!  in  (A- 
13.2) )  as 

while  Mp  and  Vp  are  obtained  from  derivatives  of  a^p,  i.e.: 


t*  4^0., 

M  =  -1 _ 2 

'  2-f  dx^ 


"  2^  dx^ 


(B-U) 


The  in  the  notation  of  Table  B-1,  the  particular  solution  to  Eq(B-l)  is 
2  "  D„  F„(X/2)-r  O,,  F„(V2)-1  f,.(V2) -f.,(X/2) 


f„(V2)-f„(A./2) 


Derivatives  of  Otp  required  to  evaluat  (B-1 5)  are  then  given  by 


At  this  point  we  wish  to  concentrate  on  evaluating  quantities  of  interest  at  (x-L)/t  =  V2  where 
the  boundary  conditions  are  applied.  We  also  assume  the  case  of  long  overlap  for  which  p,X, 
»1.  In  such  cases  all  functionaL  ratios  in  (B-1 6)  can  be  replaced  by  unity.  The  defining 
allows  Mp  and  Vp  to  be  expressed  as 

A/=-_Lci^e  ;  V  (B-19) 


-39- 


dx^  2  ’  F„(V2)-1  '^*D„F^(A/2)-1 


^21  ^12 
F„(A/2)-/-„(X/2) 


(B-16) 


dx’ 


3  -?u  /^n(^) 

‘^u  ^.M-l 


F^i%) 

F^a/2)-l 


Oh*P,y^F„(^)  -Oi, -p,)»^F„(0 
^21  ^12 
fJJT^PfjUT) 


(B-17) 


<?n  <722  ^  <7,2  fOij -»*.)'  (Pi+M2)^, 

/D,.  Y‘0^  4/  D„ 


«  ,  { 4  ,  4  m1  <722  ,  <7,2  r  (Pj  -P.)’  (P,  ^P2)’ 

Y*^..  ^023  4/  D„ 


(fl-18) 


and  comparing  (B-13)  with  (B-18)  gives  Mp  and  Vp  as  the  following  rados  with  respect  to  M„ 
and  Vo 


A/ 

P 

M. 


(B-20) 


The  values  of  these  ratios  give  a  yardstick  for  judging  how  significant  are  the  effects  of  the 
panicular  solution  on  the  peel  stresses.  NumericaJ  evaluation  indicates  that  G  and  H  are  on  the 
order  of  1  to  2.  Assuming  the  Goland  Reissner  limit  for  k  (=0.262)  applied  to  k„  (considering  the 
long  overlap  case)  and  a  value  of  1  for  G  and  H  thus  indicates  that 

=  ;  1'-0.22(eV'* 

M  ‘  V 

O  9 


The  adherend  strain  level  will  certainly  not  be  greater  than  0.01,  so  that  the  effect  of  the 
particluar  solution  on  the  peel  stresses  is  essentially  negligible. 
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